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The experimental precision that will be reached at the next generation of colliders makes it indispensable to 
improve theoretical predictions significantly. Bhabha scattering (e'*'e“ —> e'*"e“) is one of the prime processes 
calling for a better theoretical precision, in particular for non-zero electron masses. We present a first subset of 
the full two-loop calculation, namely the factorizable subset. Our calculation is based on DIANA [1]. We reduce 
tensor integrals to scalar integrals in shifted (increased) dimensions and additional powers of various propagators, 
so-called dots-on-lines. Recurrence relations [2,3] remove those dots-on-lines as well as genuine dots-on-lines 
(originating from mass renormalization) and reduce the dimension of the integrals to the generic d = 4 — 2e 
dimensions. The resulting master integrals have to be expanded to 0(e) to ensure proper treatment of all finite 
terms. 


1. Motivation 

At the next generation of linear e"*" e~ collid¬ 
ers the experimental precision of 0(1/1000) will 
exceed current standards by one order of mag¬ 
nitude. To ensure the success and strengthen 
the physics case of such a machine, a precise un¬ 
derstanding of the theoretical predictions is of 
paramount importance. In particular Bhabha 
scattering is important for luminosity monitor¬ 
ing, and hence a vital ingredient of almost every 
measurement. Recent studies [4,5] have focused 
on the calculation of the two-loop QED matrix 
element to massless Bhabha scattering, i.e. the 
photon and the electron are treated as massless 
particles. To further increase the theoretical un¬ 
derstanding of the Bhabha cross section and to 
control collinear singularities we consider the elec¬ 
tron to be a massive particle. At the energy scale 
foreseen for the next generation of e’*' e~ colliders, 
electroweak corrections cannot be ignored. While 
it is certainly a good approximation to consider 
small angle Bhabha scattering to be highly domi¬ 
nated by pure QED corrections, this does not ap¬ 
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ply to the whole phase space. Taking the electron 
mass into account in the genuine QED calculation 
will open the door to more involved multi-scale 
electroweak higher order calculations. Here we 
are presenting the first progress made towards a 
complete two-loop calculation. 

2. One-loop contributions 

At the Born level the Bhabha scattering pro¬ 
cess e“(pi)e+(p4) —> Y —> e~{-p 2 )e+{-p 3 ) is 
characterized through an s- and a t-channel pho¬ 
ton exchange diagram. One-loop radiative correc¬ 
tions constitute the two-loop calculation in the 
form of (1 loop)^ contributions to the cross sec¬ 
tion. To generate the diagrams in the form of a 
Postscript file, as well as in symbolic notation for 
further algebraic manipulation, we use DIANA, 
the Diagram ANAlyser by Tentyukov and Flei¬ 
scher [1] We choose to work in the Feynman 
gauge This has the advantage that tadpole 
contributions are energy independent and hence 
do not contribute to the physical cross section. 
Self-energy, vertex and box diagrams make their 

^For some very recent valuable features of DIANA see [6] 
as well as [7] in these proceedings. 

■^In fact, part of the calculation has been performed in a 
general gauge. 
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appearance in equal parts as s- and t-channel di¬ 
agrams. The following ten diagrams have been 
calculated. 



Figure 1. One-loop diagrams to the process 
e+e“ ^ e+e“. 


We decompose the full one-loop matrix element 
into 12 amplitudes, that is in fact into 24 ampli¬ 
tudes, counting the s- and t-channel contributions 
separately. For compactness we omit the spinors. 
The richness in structure is indeed only required 
for the box contributions. The self-energies and 
vertices contribute to only two of the structures 
below. We would like to remark at this point that 
the products of Dirac 7 matrices Mq to A4i2 can 
in principle be simplihed at the cost of artihcially 
introducing 75 , which we decided to avoid alto¬ 
gether. 

Mi = [l X 1] F’l 
M.2 = [i>i X 1] F2 

AI3 =[1 X ^2] Fa 
Mi = [i>i X ^ 2 ] F 4 
Mb = [7/i X 7;i] F 5 
Mb = [lti.il a X It,] Fq 
Mr = [it, X it,i2] Fr 
Mb = [it,iu X 7,^7^] Fb 
Mg = [ltj.li,i 4 : X 7i. 7^ ] Fg 
■Adio = [liilu X 71.74^^2] Fio 
Mii = [itj.ii,iA X 7^ 7^ i 2 ] Fii 
Mi2 = [7m7<^7p X 7p7<^7m] F12 (1) 

The entire calculation is parametrized through 
the corresponding 12 form factors for every dia¬ 
gram. Symmetries allow to further compress the 


calculation. This is in particular valuable for the 
box contributions, which are the most complex 
part of the calculation. First of all we found in¬ 
ternal symmetries relating various form factors of 
one given box diagram to other form factors of the 
same diagram: 


F2 

= F3 

Fe 

= Fr 

Fg 

= Fio 

Fi 

= 4mg Fii-1-2TOeF2 

Fi 

= 4 ITle Fio + 2 TOe Fg -f 2 TOgFg 

Fi 

= 4F8 + 2meF2 (-4 Fi 2). (2) 

The contribution in brackets is only present for 
so-called dotted-box diagrams, which will be pre- 

sented in 

the next section. 

Actually this symmetry somehow reflects the non- 
minimal choice of amplitude-basis we made. Fur¬ 
ther conventional crossing symmetry can be ap¬ 
plied to relate s- and t- channel diagrams. More¬ 
over we found one additional symmetry relating 
crossed box (cb) diagrams to direct box (db) di- 

agrams in 

the same channel: 

F^'^ = 

Ff- 2me Fg^*’ + 2 dFt- 4 :vie F^ 

cr 

II 

-Ft -‘ 2 . d Ft 

cr 

II 

Ft -{ 4 - 2 d) F{^i‘’ 

Ft = 

-Ft + (4 - 6 d)Ft 

Ft = 

-Ft -4 Ft 

Ft = 

-Ft 

cr 

II 

Ft 

II 

1 

cr 

II 

Ft, (3) 


and interchanging Ft'°{s,t,u) Ff'°{s,u,t) suc¬ 

cessively. The space-time dimension d reflects our 
choice of dimensional regularization. It is suffi¬ 
cient to calculate one box diagram and obtain the 
form factors of the remaining three box diagrams 
via those symmetry relations. 

To reduce the various one-loop integrals we rely 
upon two different techniques. First we follow the 
prescription of Passarino and Veltman [ 8 ] with 
the subtle difference from the original work that 
the tensor decomposition is performed with re¬ 
spect to inner momenta. 
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Favourably we use the method first proposed 
by Davydychev [9] to reduce one-loop tensor inte¬ 
grals into scalar integrals and changing the space- 
time dimension d, and later extended to multi¬ 
loop tensor integrals by Tarasov [2]. Let 
be a one-loop n-point vector (tensor) integral, k 
the loop-momentum, Cr = {k — ~ nT-r the 

propagator corresponding to the r-th line of the 
diagram appearing in the i/^-th power, with the 
special case c„ = Those integrals can 

be written as 




TflU 



n—1 


- ^ pf (1 + S,,) i g>^’' /[f+l, 

*,i=i 


where [c?-|-] is an operator shifting the space-time 
dimension by two units and 

= / n Ur + Sri+Srj- 

J r=l Cr 

is the original scalar integral with one additional 
power in the i-th and j-th propagators. Hav¬ 
ing reduced the tensor integrals to scalar inte¬ 
grals the generic space-time dimension d needs 
to be re-established. For this we use the recur¬ 
rence relations first proposed in [2], which are 
complementary to those obtained via integration 
by parts [10,11], and later simplified and ex¬ 
tended to zero Gram determinants in [3]. Let 
Yij = — (pj — pjY + mf + m'j and the Cayley de¬ 
terminant 



0 

1 

1 . 

1 


1 


^2 . 

.. Yin 

On ^ 

1 

Yi2 

^22 

.. Y2n 


1 

Yin 

Y2n ■ 

Y 

• • J nn 


then the so-called signed minors (^^ are de- 
fined as the rows ji,j 2 ,... and columns fi,Z 2 ,... 
erased from the Cayley determinant ()n-^ Making 

^Note here the additional overall sign resulting from 
(_l)Ji+i2+...+n+i2+.., 


successive use of the following three recurrence re¬ 
lations leads to scalar master integrals Ao,Bo, Cq 
and Dq in d dimensions: 


On 


+ j{d-\-2) 
■^n 



rid) 

^ TJ. 




+ rid) _ 


- E 

i,k 


Pi \ 


Ad) 


1 -f - d 



n 


E 



{Vk - 1 ) 


Ad) 

■^n 





the operators i^, raising the power of the cor¬ 
responding propagator by one unit , while k~ re¬ 
duces the power of the fc-th propagator by one 
unit. For zero Gram determinants ()„ = 0, simi¬ 
lar recurrence relations are needed and are avail¬ 
able in [3]. Effectively a zero Gram determinant 
reflects the kinematical boundaries of phase space 
where a given n-point function can be expressed 
through scalar integrals of lower rank. Typical 
examples of such simplifications are 




{d-2) 1 Ao{ml) 


(d — 4) s — Ami 


+2 


(d-3) _ 1_ 

(d — 4) s — 4? 


■Bo{s;ml,ml) 


Bo{ml;ml,0) 


1 (d-2 ) 24o(rn^) 

2 (d — 3) TOg 


3. Factorizable two-loop contributions 

Two-loop contributions to the Bhabha matrix 
element can be classified according to the topo¬ 
logical properties of the corresponding diagrams. 
In this section we focus on the contributions orig¬ 
inating from the renormalization of the electron 
mass. To accommodate the renormalization con¬ 
tributions, every internal electron propagator in 
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Fig. 1 has to be replaced by 

1 ^ A _|_ ^TOgl^TOe 

PI — {me + 5meY Pe ~ V Pe ~ 

with 


( 10 ) 


5me = —(1 — e) 

-2 e'^ me Boipl;ml,0) 


me 

„ 2 . ^2 


the electron mass counter-term. From Eq. (10) we 
understand that we have to calculate all possible 
permutations of the diagrams in Fig. 1 with one 
of the loop-internal electron propagators squared. 
Now, since the mass renormalization term 5m^ 
contains 1/e UV poles, the one-loop master inte¬ 
grals Aq,Bq,Cq and Dq are needed in 0(e) to 
ensure that all finite terms are taken into ac¬ 
count properly. The e-expansion of Aq is pretty 
straightforward. The hypergeometric presenta¬ 
tion of 1-loop scalar master integrals in arbitrary 
dimensions d was given in [12,7]. From these re¬ 
sults the e-expansion can be obtained and will 
be presented in [13]. The electron propagator 
appearing squared does not really impose fur¬ 
ther complications in the recurrence relations ap¬ 
proach. For this part of the calculation we refrain 
from further pursuing the independent calcula¬ 
tion based upon the Passarino-Veltman reduction 
method. 


4. Results 

Given the scope of this article we cannot 
present the full analytical result here To give 
an idea of the structure we present exclusively the 
form factors F/ for the various s-channel topolo¬ 
gies in the Appendix. We use the following ab¬ 
breviations 


^0 

= Ao{ml), 


Bs 

= Bo{s;ml,ml), 

Bt = Bo{t;ml,ml) , 

Co 

= Co{ml,ml,s;0, 

ml,0) , 

Do 

= Do{me,me,me, 

ml,t,s;ml,0,ml,0), 


= l/{l-Aml/s), 

= 1/(1 _ Aml/t ), 


drUp t 

u 


= , rt = 

, ?■„ = -. (11) 


s s 

s 

4 For 

a complete result see 

http://acat02.sinp.msu.ru 
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5. Summary 

In this article we have presented the complete 
one-loop originating contribution to the two-loop 
Bhabha cross section. We decomposed the matrix 
element into 12 amplitudes and presented the re¬ 
sult in the form of the corresponding form factors. 
Mass renormalization lead to dotted diagrams, 
which we reduced to scalar integrals by making 
use of recurrence relations. Two-loop factorizable 
form factors were presented. 
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A. Appendix 

All form factors have the overall factor 
e^/(47r)‘^/^. Form factors corresponding to dot- 
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ted diagrams have the additional common factor 
Both box contributions have the addi¬ 
tional common factor 

^dotself ^ Ao\{d-2)x^z^ 

-Bs-z^[l-{d-3)x^] (A.2) 

s 




= 2 { -A, 


:[2 


d- 3 


2 2i 

— X Z \ 


\mi 

+ Bs^[x^{l -I- -I- d - 4] 


(d — 4)s 


[1 + ^^] 


K 


'dot vertex 


= 2 


-f 6- 


d-3 d-5 


(A.3) 
1 

:) 


-(d - 4){2x^ -l + {d-4:))+7- 4x^] 


- Bs-x^z^[{d - 4:){2x^ + {d- 4)) -f 2x^] 


An 


+ B, 


(d-4), 


2 2 
-X Z 


(A.4) 


-rf + nrux'^) + z'^irtXu - n + irl - TuV^ 

-z'^rt)) - (2rt(l + rt + 2r„ -h rtr^x^) 

+2z^{l + ru)y^))Bt — {rt{l +rt + Xu) 
s 

+z‘^{l + ru)y^) - Bs — {rtru{l + nx^)) 
s 


A, 


0 


+ Bt 


2n 


-{2rt{rt +rl 


s{d — 4) 

+ {2 + rt)ru) + z^{rt{l - r„ - 3r0 


-f(2 -h 3r„)y2 + z^n)) - 

(d-3) 

x(2(rt - r? (1 - rl)ru) - z'^{2 + r* -f r„ 
-3rtr„ - 5rj - z'^{l - r„ - 4rt + z^))) 
-2{rt + rf + (1 + 2rt)ru - z‘^rt{l - r„a;^))) 


FDortsi ^ A rt - + z^{rt + 3r^ - 2r^ 

(d-3) 

-1-2(1 - r'^)ru - z^(2(l + n) - br^ + (1 - 3rt)r„ 
-z^(l - r„ - 4rt -h z^)))/2) -h rt(l -f rt){l + rt 
+2ru) - z^(4rt(l -hrt)-(r„ -4rt)r„ - 2z\)l2) 

(A.5) 


J^dotdbox ^ Ao^{{d-A){2r^{r‘^+irt + 2 

m^s 

+ru{2 + 2rt + 2r„ - x'^rt)) - z^{rti3rf + brt 
+2rtru - y^{l + 2ru)ru) - z^{rt{rt - y^) 

-(1 +?’«)y‘‘))) + j^j^^{2r‘l{r‘l+ 2,rt + 2 

+ru{2. -f 2r„ -F rt)) - z^rti^rl + 3rt + rtVu 
-y‘^2{l + ru+ rl) - z^rt)) + 2r^{3r^ + 9rt 
+5rtr„ -h 6(1 + r„ -h rl) - 2x^rtr„) -h z^(rt(l 
+ru - 22rt - 14rt - 6rtr„ +y'^{3 + 17r„ 

+ 12r2))/2 + z^irt{8rt - 2 - 4y^{l - r„)) 

-4y\2 + r„))/4)) - Bt-(2(d - 4)(2rt(rt 
s 

+y‘^il + Tu)) + z^y‘^{rt + y^(l + r„))) - rt(l 
-|-4rt + (1 + 4rt)r„ - 5y^(l - r„)) -|- z^(rt(l + 

2rt) - 2j/2(rtr„ - y'^))) - Bs-{d - 3){8r'^{l + rt 

s 

+ru)ru + z^rt(l -|- 3rt -|- ru(l - 4rt) -|- 2x^rt(l + 
2rt)r„ -I- 2y‘^{l + 2r„ -|- 2rl) - z^(l -|- 2x^rtr„))) 

+ (- 2 F Ft] —-^(rt(l + Xu + 4rt(l + Xu 

\me^ /yd — 4)s 

+ft) - 2/^(1 - 9r„)) - z^(rt(l + 2xt - 2y^(l -|- 
r„)) - 2y‘‘r„)) -I- CoXt{{d - 4)(4rt(l + rt -I- r„)ru 
-|-z^((l + r„ - rt(l + 10r„ -|- 8r^ -I- rt(8 -I- 8xu + 
4rt)) + 4a;^rt ru -h 2y'^{l + 2r„ -h 2r^))/2 - z^(l 
-4rt(l + Xu) - 6xt + 2/^(2 -f 4xu + 4r^) -f 2z^rt) 
/2)) - z^(rt(r„ -h 2rt(l -h rt) - 2x^rtr„) - z^(l 
+Xu + 3rt -I- 6rt - z^(l + 2rt))/2)) -|- DoXtz'^s{{d 
-4)(rt(l + 2(rt -h (1 + rt + r„)r„) -h r^) - z^(rt 
(2 -h 2r„ -h 3rt) - 2/^(1 -h 2xu + 2x1) “ z'^rt)l‘2.) 
—Xt{2, 2xt + 2rtXu + J'u + 8/2xu) + z^(rt(3 -I- 
2r„)-2/"(l + 2r„ + 2r2))/2) (A.6) 



